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ABSTRACT 



Due to the recent dramatic technological advances, infrared interferometry can now be applied to new classes 
of objects, resulting in exciting new science prospects, for instance, in the area of high-mass star formation. 
Although extensively studied at various wavelengths, the process through which massive stars form is still only 
poorly understood. For instance, it has been proposed that massive stars might form like low-mass stars by mass 
O ' accretion through a circumstellar disk/envelope, or otherwise by coalescence in a dense stellar cluster. Therefore, 
-£3 , clear observational evidence, such as the detection of disks around high-mass young stellar objects (YSOs), is 
^ ' urgently needed in order to unambiguously identify the formation mode of the most massive stars. 

After discussing the technological challenges which result from the special properties of these objects, we 
, present first near-infrared interferometric observations, which we obtained on the massive YSO IRAS 13481-6124 
' using VLTI/ AMBER infrared long-baseline interferometry and NTT speckle interferometry. From our extensive 
, data set, we reconstruct a mo del- independent aperture synthesis image which shows an elongated structure with 
■ a size of ~ 13 x 19 AU, consistent with a disk seen under an inclination of ~ 45°. The measured wavelength- 
dependent visibilities and closure phases allow us to derive the radial disk temperature gradient and to detect a 
dust-free region inside of 9.5 AU from the star, revealing qualitative and quantitative similarities with the disks 
observed in low-mass star formation. In complementary mid-infrared Spitzer and sub-millimeter APEX imaging 
observations we detect two bow shocks and a molecular outflow, which are oriented perpendicular to the disk 
plane and indicate the presence of a bipolar outflow emanating from the inner regions of the system. 
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1. INTRODUCTION 

In recent years, the field of infrared interferometry has experienced dramatic technological progress, providing 
substantial improvements, for instance, in terms of the available baseline lengths, spectral resolution, and spec- 
tral coverage. Besides quantitative improvements, the instrumentational and infrastructure advances have also 
considerably extended the accessible range of astrophysical targets, enabling infrared interferometry to address 
key scientific areas such as high-mass star formation. 
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Although scarce in number, massive stars (> 10 Mq) are of fundamental importance for astrophysics. Already 
shortly after their birth, they start to disperse their natal molecular clouds with strong stellar outflows and 
photoevaporate the protoplanetary disks around nearby low-mass stars. After a short hydrogen-burning phase, 
the most massive stars end their evolution as supernovae, enriching the interstellar medium with heavy elements 
and triggering the next generation of star formation through supernova shock waves. In spite of their significance, 
it is still poorly understood how high-mass stars form. 

Early spherically symmetric calculations suggested that the accretion scenario, which appears well established 
for stars with masses below 10 Mq, might not work for massive protostars due to the strong radiation pressure^^. 
Therefore, it was proposed that high-mass stars might not form by accretion, but instead by stellar merging^. 
Although recent theoretical work showed that the radiation pressure limit might be overcome when considering 
more complex infall geometries than spherically symmetric onea^H clear observational evidence, such as the 
detection of compact dusty disks around massive young stellar objects, is needed to unambiguously identify the 
formation mode of the most massive stars. 

Intcrferometric observations at radio wavelengths already provided fascinating insight into the distribution 
and kinematics of the circumstellar material, although the existing observations are still strongly limited by the 
achievable angular resolution of thousands to hundreds of astronomical units (AUs). For instance, millimeter 
observations^ with an angular resolution of ~ 0.8" revealed a large (~ 10 4 AU) rotating toroidal structure and gas 
infall motion around the YSO G24.78+0.08. Recently, mid-infrared long-baseline interferometric observations^"^ 
probed the environment around two massive YSOs on scales of a few tens of AUs and could resolve the inner 
regions of an envelope around these stars. In order to separate the contributions from a circumstellar disk and to 
avoid confusion with the surrounding rotating circumstellar envelope, observations in the near-infrared around 
~ 2 /im are very promising, since this wavelength regime is sensitive to the thermal emission of hot (> 1000 K) 
gas and dust, probing temperatures where dust is expected to sublimate (T ~ 1500. ..2000 K) and to form a 
dust-free inner hole. By forming telescope baselines up to hundreds of meters, infrared intcrferometric imaging 
can now provide the milli-arcsecond (mas) angular resolution that is required to directly resolve these inner 
disk regions. Compared to earlier studies on high-mass YSOs, which used conventional imaging techniques at 
10 m class telescopes equipped with adaptive optics systems, the gain in resolving power is at least one order 
of magnitude, while the gain compared to state-of-the-art (sub-) millimeter disk studies is about two orders of 
magnitudes. 

However, the technological challenges for infrared interferometry on young high-mass stars are high. First, 
massive YSOs are still deeply embedded in their natal clouds, which results in an extremely steep drop of the 
spectral energy distribution (SED) between mid-infrared and near-infrared wavelengths. Accordingly, most mas- 
sive YSOs are not associated with a visual counterpart. Therefore, it is essential to use off-axis telescope guiding, 
as offered by the STRAP units at VLTI 1.8 m auxiliary telescopes (ATs), which use Avalanche Photodiodes to 
provide a tip-tilt correction for guide stars brighter than V = 13.5. The VLTI 8.2 m unit telescopes (UTs) are 
equipped with MACAO adaptive optic units operating on natural guide stars down to V = 17. These sensitivity 
limits and the requirement that the guide star has to be located within the isoplanatic patch (~1'), makes it 
difficult to find suitable guide stars for many high-mass YSOs, in particular since the number of potential off-axis 
guide stars is also strongly reduced by the high visual extinction in high-mass star forming regions. 

Recently, we have obtained first near-infrared interferometric observations of a young high-mass staiP, namely 
IRAS 13481-6124. This star has a nearby foreground star (V — 12.6) at a separation of ~ 17", which can be 
used for telescope guiding even with the ATs. Furthermore, IRAS 13481-6124 is located at a very favourable 
celestial location (DEC=— 61°) which results in rather circular uu-tracks, providing optimal prerequisits for 
earth-rotation aperture-synthesis imaging. Earlier studies have estimated that IRAS 13481-6124 is located at a 
distance of about 3.5 kpcP^l and harbours a central object with a mass of ~ 20 Mgpl, which is embedded in a 
cloud with a total mass of ~ 1500 MgP. 
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Figure 1. Left: un-plane coverage obtained with our VLTI/ AMBER and NTT/speckle observations. Right: Visibilities 
recorded with our VLTI/ AMBER and NTT/speckle observations towards different position angles, revealing a strong 
position angle-dependence of the visibility function. 

2. OBSERVATIONS 

We observed IRAS 13481-6124 during seven nights in 2008 and 2009 using the near-infrared three-telescope 
beam combination instrument AMBER 13 at ESO's Very Large Telescope Interferometer (VLTI). Our spectro- 
interferometric observations cover the near-infrared if-band between ~ 1.95 to 2.55 /im with a spectral resolution 
A/AA = 35. Using the VLTI/ATs on three different array configurations (E0-G0-H0, D0-H0-G1, A0-D0-H0) 
and at various hour angles, we could achieve a relatively good uu-coverage (Fig. [TJ left), including baseline 
lengths between 12 and 85 m. In order to cover spatially extended structures, we complemented the VLTI 
long-baseline interferometric observations with bispectrum speckle interferometric observations^^, yielding 
precise visibility information for baselines < 3.5m. The speckle if-band data was recorded with ESO's New 
Technology Telescope (NTT) and our visitor camera, employing a Rockwell HAWAII detector. Following the 
standard observing procedure, both the VLTI and speckle observations on IRAS 13481-6124 were interlayed with 
calibrator observations in order to derive absolute calibrated visibilities and closure phases, which we then used 
for our image reconstruction and model fitting approach. 

Besides these high-resolution if -band continuum observations, which directly trace the thermal dust emission 
and scattered light in the circumstellar environment around IRAS 13481-6124, we searched for signatures of mass 
outflow on larger spatial scales. For this, we obtained archival Spitzer/IRAC images which were recorded in the 
course of the GLIMPSE survejE^ and which cover wavelengths around 3.6, 4.5, 5.8, and 8.0 /im. In the IRAC 
4.5 /xm image (codes as green color in the color composite shown in Fig.[2K) we detect two bow shock structures 
which are oriented along position angle 31 ± 6° and likely indicate H2 shock-tracer line emission^. In order to 
trace a potential molecular outflow, we obtained a 12 CO (3-2, 346 GHz) map using the SHFI heterodyne receiver 
at the APEX 12 m telescope located on the Chajnantor plateau in the Chilean Andes. The In the obtained 
CO channel maps, we detect a bipolar outflow oriented along position angle 26 ± 9°, with the approaching 
(blue-shifted) lobe southwest of IRAS 13481-6124 (Fig.[2p). 

3. APERTURE SYNTHESIS IMAGING 

Given that our observations provide a good uw-coverage and include closure phase information, we could use our 
data to reconstruct a model-independent aperture synthesis image. For this purpose, we employed the Building 
Block Mapping (BBM) image reconstruction algorithm 18 . This algorithm has already been employed successfully 
in various earlier projects to reconstruct interferometric images from simulate d 19 !! 20 ^ and long-baseline interfero- 
metric data^i" 2 ^. The reconstruction was performed on a 256 x 256 grid using a pixel scale of 0.85 mas/pixel and 
convergence was reached after ~ 12 000 iteration steps. In order to maximize the uu-coverage for imaging, we 
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Figure 2. Zoom in on IRAS 13481-6124, covering spatial scales over more than five orders of magnitudes: (A) In 
Spitzer/IRAC images, we detect two bow-shock structures, indicating a collimated outflow. (B) The outflow is also 
detected in molecular line emission using APEX/SHFI on scales of a few 10 000 AU (the contours show blue-shifted veloc- 
ities of -10 km s~ 1 and redshifted velocities of +5 km s~ 1 , where the latter are indicated with dashed lines) . (C) Applying 
our image reconstruction algorithm to the VLTI/AMBER data, we reconstructed a model-independent aperture synthesis 
image of IRAS 13481-6124. The image reveals an elongated structure, which is oriented perpendicular to the outflow 
direction (contours decrease from peak intensity by factors of y/2). (D) Elongation of the compact emission component, 
as determined with our temperature gradient disk model. 



combined data from all spectral channels in order to reconstruct one monochromatic image. The obtained image 
was convolved with a Gaussian of FWHM 2.4 mas, corresponding to a resolution of X/2B, where A denotes the 
observing wavelength and B the length of the longest employed baseline. 

The reconstructed image (Fig.[2p) has a reduced goodness-of-fit value of 1.5 and 1.0 for the squared visibilities 
and the closure phases, respectively. The residuals are well distributed, i.e. they closely follow a normal distri- 
bution. The image clearly reveals an elongated structure which is oriented in north-western direction. Besides 
this compact component, some flux appears distributed over the whole image. Measuring the flux in different 
apertures, we find that 80% of the total image flux is contained in a 7.0 mas aperture, while the remaining flux 
elements are spread rather homogeneously over the image. 

4. MODELING 

In order to further characterize the elongated structure detected with aperture synthesis imaging, we fitted 
geometric and detailed physical models to our data, allowing us to deduce object information on smaller scales 
than the diffraction-limited resolution of our reconstructed image. 

4.1 Geometric models 

Our measurements show that the visibility function drops rapidly at baseline lengths < 3 m, but stays then nearly 
constant up to baseline lengths of about 12 m, followed by a uniform, nearly linear decline at baseline lengths 
up to 85 m (Fig. [TJ right), providing clear evidence for the presence of at least two spatial components. We use 
a Gaussian model to estimate the characteristic size of the two components and find that the first, extended 
component has a FWHM diameter of about 108 mas and contributes ~ 15% of the total flux. Based on the 
finding of scattered light envelopes around other massive stars, we suggest that this extended, centro-symmetric 
component, which also appears in our image as homogenously distributed flux elements, is the scattered light 
contribution of a natal envelope. The second, compact component accounts for 85% of the total if -band flux 
and has a size of 5.4 x 3.8 mas, consistent with a disk seen under an inclination angle of 45°. Furthermore, we 
find that the size of the compact component increases significantly towards longer wavelengths, as shown by the 
measured equivalent widths in Fig. [3] (top). This effect, which was found towards many low- and intermediate- 
mass YSO disk a * * \ is likely indicating a temperature gradient in the circumstellar material, where hotter 
material (radiating more effectively at shorter wavelengths) is located closer to the star. 

Motivated by these indications for an internal temperature structure, we fitted analytical disk models with a 
radial temperature power-law, T(r) = Xi n (r/ri n ) -9 , to our data. Assuming reasonable values for the inner disk 
temperature, i.e. T- ln = 1500 to 2000 K, we find that this model can reproduce our interferometric data with an 
inner disk radius r- ln of 3 mas (9.5 AU), and a temperature power law index q of 0.4, which is consistent with the 
theoretical temperature gradient of flared irradiated disks^l (q = 0.43). Compared to other geometric models 
(e.g. uniform disk model with \ 2 = 4.48 or Gaussian model with \ 2 — 2.56), the temperature gradient disk model 
provides a significantly better representation of our data (\ 2 — 1.40). Intriguingly, the derived inner disk radius 
of 9.5 AU agrees with the expected location where dust in an irradiated circumstellar disk would sublimate, i.e. 

6.2 AU to 10.9 AU assuming grey dust and dust sublimation temperatures between 2000 and 1500 K. Therefore, 
IRAS 13481-6124 is following the size-luminosity relation^, which is well established for low- to intermediate 
mass YSOs, suggesting that the near-infrared emission mainly traces material at the dust sublimation radius, 
similar as in the disks around T Tauri and Herbig Ae/Be stars. 

The derived disk orientation (120°, Fig. 03) is perpendicular to the determined outflow axis (26 ± 9° and 
31 ±6°), suggesting that the AU-scale disk resolved by VLTI/ AMBER is indeed the driving engine of the detected 
collimated outflow. 

4.2 Radiative transfer modeling 

In order to confirm that the proposed scenario, including a high-mass young stellar object, a compact disk, and 
an extended circumstellar envelope, can fit not only our interferometric observations but is also consistent with 
the measured SED, we performed detailed radiative transfer simulations. To obtain an initial parameter set 
for our modeling procedure, we used the pre-computed grid of 2-D radiative transfer models by Robitaille et 



dRA [AU] 
10 5 5 10 



O 

LU 
Q 

"O 



2.5 iim 
2.3 urn 
2.1 urn 




Gaussian profile 



10 




5 










< 





O 




LU 




Q 




"O 


5 




10 





2 2 

dRA [mas] 




o 

1.0 

0.8 

0.6 

0.4 

0.2 


1.0 

0.8 

0.6 

0.4 

0.2 







k, 100°. 

V 




120° . 


k K 130°; 


"K 140 °- 







5 



10 15 20 25 30 35 5 10 15 20 25 30 35 5 10 15 20 25 30 35 5 10 15 20 25 30 35 
Spatial freq. [10 6 X] Spatial freq. [10 6 X] Spatial freq. [10 6 X] Spatial freq. [10 6 

Figure 3. Elongation of the near-infrared emitting region (top), as determined by fitting Gaussian models to independent 
subsets of our visibility data (bottom), each covering a position angle range of 10°. In order to investigate whether the 
object morphology changes with wavelength, we have divided our data in three wavelength bins (blue: A < 2.2; green: 
2.2 < A < 2.4; red: A > 2.4) and then adjusted the Gaussian FHWM diameter. Besides the compact component, this 
model also includes a Gaussian with FWHM diameter 108 mas, representing the extended envelope component detected 
with our NTT/speckle observations. 
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Figure 4. Top: Model image (left, A = 2.25 (im) and SED (right) corresponding to our 2-D radiative transfer model. 
Bottom: Comparison of the measured and the model visibilities for the wavelength range covered by our AMBER obser- 
vations (1.95 to 2.55 fim, grey dashed lines). The black lines show the visibility function for the central wavelength of 
2.25 (j,m. 



alPSl. The model parameters are then put into an adapted version of the radiative code by Whitney et alP3 
in order to compute model SEDs and synthetic images, from which we compute visibilities and closure phases 
for comparison with our interferometric observations. After some parameter adjustments we find a parameter 
set which reproduces the SED and the measured visibilities and closure phases reasonably well (Fig. 0]). The 
model includes a central star with a mass of ~ 18 Mq, a circumstellar disk with a curved rim and an inner dust 
truncation radius at 6.2 AU, and an extended envelope with bipolar outflow cavities. For a detailed description 
of the model components and the best-fit parameters, we refer the reader to our recent science papei^. 

5. CONCLUSIONS 

Probing the dust distribution around a young high-mass star on the smallest physical scales yet, our study 
provides compelling evidence for the accretion-disk hypothesis in high-mass star formation and permits a first 
characterization of the geometry of the innermost disk regions. Furthermore, our study demonstrates the benefits 
of a multi-wavelength (2-870 /im) and multi-scale imaging approach, which allows us to relate the detected 
AU-size compact disk with the outflow signatures detected on parsec scales, yielding a global picture of the 
environment around this massive young star. 

With the available near-infrared interferometric instrumentation, the number of accessible high-mass YSO 
targets is rather limited, in particular due to the current sensitivity limits and the need for off-axis telescope 
guiding. However, the number of accessible targets might increase considerably with the upcoming generation 
of infrared interferometric instruments. For instance, the implementation of infrared wavefront sensing units 
would enable interferometric observations on a large number of targets without suitable visual guide stars. For 
VLTI, this implementation of near-infrared wavefront sensors (H- or JVband) is planned in the course of the 
GRAVITY projeclP. 

Besides these improvements for telescope guiding and beam injection, further progress could be achieved 
by observing in an optimal wavelength regime. Near-infrared (H- and iVband) interferometric observations 
provide an excellent angular resolution and already enable phase-closure imaging observations, but are restricted 
to a small sample of targets. In the mid-infrared iV-band (8-13 /im), the brightness limitations are significantly 
relaxed, while at these wavelengths the disk contributions get more contaminated with emission from the envelope 
and the angular resolution becomes insufficient to study the inner disk truncation region. Therefore, it is 
likely that the optimal observing wavelength for many young high-mass stars is between the K- and A-band, 
for instance in the L- or Af-band. Two-telescope L-band interferometry was recently demonstrated at the 
Keck Interferometei)^ and will also become available with the VLTI second-generation instrument MATISSE^], 
enabling four-telescope L- or Af-band interferometric imaging on a large sample of massive YSOs. 
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